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ABSTRACT 



In 1998, the Super-Kamiokande announced evidence for the observation of noii- 
trino oscillations based on measurements of the atmospheric neutrino fluxB'. 
This paper presents the updated results for fully and partially-contained events 
with 736 days of data, as well as upward-going muon results and a global anal- 
ysis. Preliminary interpretations of the results in terms of various two-flavor 
oscillation hypotheses are presented. 



1. Introduction 

1.1. Atmospheric Neutrinos 

Atmospheric neutrinos are produced by collisions of cosmic rays with atoms in the 
Earth's upper atmosphere. These collisions produce showers which include neutrinos 
with energies in the 0.1 GeV to 100 GeV range. Although calculations of the absolute 
flux of atmospheric neutrinos have fairly large uncertainties (~ 20%), the ratio of 
muon to electron flavor content is much better known (to ~ 5%), since it arises from 
the decay chain ir — > /iz/ M ,/i — > ev^ v e . The muon to electron flavor ratio is expected 
to be about 2 (one v e from muon decay for every two u^), gradually rising at energies 
above about 5 GeV as the muons from ir decay have enough energy to survive to the 
surface of the earth and range out before decaying. 

Because the attenuation of even high energy neutrinos in the earth is so small, 
neutrinos produced all around the atmosphere can be detected. Their pathlengths 
range from about 15 km for neutrinos produced directly overhead, to 13000 km for 
neutrinos produced directly below. The pathlength of the neutrino can be inferred 
from the direction of the Cherenkov cone of its produced lepton. Therefore, by 
measuring the zenith angle, energy and flavor of atmospheric neutrinos, we can test 
the neutrino oscillation hypothesis, which gives the probability of flavor conversion: 




(1) 



1.2. The Super- Kamiokande Detector 



The Super-Kamiokande (Super-K) experiment is a large water Cherenkov detector 
located in Mozumi, Japan. Its total mass of ultra-pure water is 50 kton, divided into 
two concentric cylinders: an inner volume with its inside surface covered by 11146 
inward-looking 50 cm photomultiplier tubes, and an outer volume serving as entering 
particle shield and veto with 1885 outward-looking 20 cm phototubes. The fiducial 
mass of the inner volume (2 m away from the walls) is 22.5 kton. 

2. Atmospheric Neutrinos with Super-K 

High energy neutrinos produced in the atmosphere are detected via their interac- 
tions with nuclei in the water. The easiest v interaction to tag and study in a water 
Cherenkov detector is the charged current (CC) quasi-elastic interaction: 

ui + N — I + N', (2) 

where the flavor I of the outgoing lepton tags the flavor of the incoming neutrino. The 
flavor of this lepton is determined using the Cherenkov ring: events are identified as 
"//-like" (non-showering) or "e-like" (showering) based on a likelihood analysis of the 
Cherenkov light around the cone. At higher energies (more than a few GeV), more 
complicated single and multi-pion products from both charged current and neutral 
current (NC) neutrino interactions with nuclei predominate. 

Atmospheric neutrinos are divided into three classes for Super-K: 

1. Fully-contained events (FC): events for which the products of the neutrino 
interaction are completely contained within the inner detector volume. These 
are further subdivided into "sub-GeV" (with visible energy <1.33 GeV) and 
"multi-GeV" (>1.33 GeV) events. For the current analysis, single ring (mostly 
quasi-elastic charged current) interactions are selected. According to Monte 
Carlo (MC), the FC e-like sample contains about 86% CC v e , 4% CC and 
10% NC interactions; the FC //-like sample contains about 96% CC u^, 4% CC v e 
and <1% NC interactions. 

2. Partially-contained events (PC): events for which the produced muon exits 
the inner detector volume (only muons are penetrating enough). Such events 
comprise a 98% pure sample of CC interactions, according to MC. There is 
no single ring selection for partially contained events. The average energy of a 
neutrino producing a PC event in Super-K is about 15 GeV. 

3. Upward-going muons: upward-going muons result from neutrinos which in- 
teract in the rock below the detector to produce a penetrating muon. As for PC 



events, we assume the neutrino parent flavor is u^, since only muon flavor neu- 
trinos can produce penetrating leptons. These muons may either traverse the 
entire detector or stop inside the volume. Through-going upward muons have 
average parent neutrino energies of about 100 GeV; stopping upward muons 
have average parent neutrino energies of about 10 GeV. (Of course, there are 
also downward-going neutrino-induced muons, but we cannot distinguish these 
from the large background of entering downward-going cosmic ray muons). 

Details of the event selection and reconstruction (energy, vertex position, direc- 
tion, and number of rings), and of the MC are described elsewhere 

3. Recent Results 



3.1. What's New 

The new preliminary data sample, approved in December 1998, is described in 
Table 0. 

Table 1. Summary of June 1998 and December 1998 data and MC samples. 



DATA 




June 1998 


December 1998 


Analyzed ivetime: 






Fully-contained 


535 days 


736 days 


Partially-contained 


535 days 


685 days 


Upward-going muons 


516 days 


516 days 


MONTE CARLO 




June 1998 


December 1998 


Generated livetime 


10 years 


20 years 


High energy cross-sections 


CCFR 


GRV94 PDF& 


Solar epoch 


solar cycle average 


solar cycle minimum 



This new sample incorporates additional livetime for FC and PC events. There 
were some minor improvements in reconstruction (fitting, ring counting, etc.). The 
MC, however, is a completely new sample. The most important differences are in- 
dicated in the table. The net effect of the change is an increase in the number and 
the fraction of multi-ring events at higher energies. Twenty years of livetime were 
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Fig. 1. The east-west effect: azimuthal distribution for low-energy near-horizon e-like and /i-like 
data, along with MC predictions based on Honda (solid line) and Bartol (dashed line) fluxes. <j> = 
corresponds to particles going to the north and <j> = tt corresponds to particles going to the south. 
The expected deficit of neutrinos from the east due to the geomagnetic cutoff of the charged primaries 
is observed. 



generated (previously there had been only 10 years). 



3.2. Fully and Partially- contained events 



3.2.1. East- West Effect 



The east- west effect results from the action of the geomagnetic field on the charged 
cosmic ray primary particles. Protons at the horizon with energies less than about 
50 GeV/c coming from the east are suppressed relative to those coming from other 
directions, where the momentum cutoff is smaller (~10 GeV/c). To check the relia- 
bility of both our experimental results and the atmospheric neutrino flux calculations, 
we have selected a sample of neutrinos for which the east- west effect is expected to be 
important: for which |cos 6 1 ! < 0.5 and momentum is between 400 and 3000 MeV/c. 

Figure [l] shows the azimuthal distribution for both e-like and /z-like events, along 
with two flux predictions!'!^. The expected deficit from the east is observed. This 
analysis has been described in detail elsewhereS-* . 



3.2.2. Ratio of Ratios 



The measured flavor ratio is conventionally compared to expectation as the "ratio 



of ratios" , R, defined as 

„ (NjN e ) DATA 
{NjN e ) MC ' 

where and N e are numbers of /i-like and e-like events for data and MC 
numbers (Honda flux!-*) are given in Table [^, and the new values of R are 

R FC = 0.67 ± 0.02(stat) ± 0.05(sys) (4) 

for sub-GeV FC events, and 

R FC+PC = 0.66 ± 0.04(stat) ± 0.08(sys) (5) 

for multi-GeV FC and PC events. If there are no neutrino oscillations, R is expected 
to be 1. However, R is significantly smaller than 1 for both sub-GeV and multi- 
GeV data. A significantly low value of R has been measured by several experiments 
in the past00\ and recently by the Soudan collaboration with an iron-calorimeter 
detector©. 

Tabic 2. Numbers of events for FC and PC data and Monte Carlo. 



(3) 
. These 



SUB-GEV 




Data 


MC 


Single ring: 






e-like 


1607 


1510.5 


/i-like 


1617 


2277.8 


Multiple ring 


1271 


1614.3 


Total 


4495 


5402.5 



MULTI-GEV 




Data 


MC 


Single ring: 






e-like 


386 


357.4 


/i-like 


301 


415.9 


Multiple ring 


737 


925.7 


Total 


1424 


1698.9 


PC (/i-like) 


374 


528.7 



3.2.3. Angular Distributions 

To study the pathlength dependence of the disappearance of atmospheric muon 
neutrinos, we can study the flux as a function of zenith angle: the Cherenkov ring 
direction is correlated with the neutrino direction. Angular distributions for FC and 
PC events are shown in Figure |2|; cos 8 = — 1 corresponds to upward-going neutrinos 
and cos 9 = +1 corresponds to downward-going. The oscillation prediction for the 
best fit parameters (see the next section) for — > v T oscillation is superimposed. 

For multi-GeV events, Figure ^| shows the angular distributions in 10 bins, showing 
that the oscillation prediction matches the data through the rapidly changing region 
near the horizon. 

The zenith angle asymmetry A, defined as A = {jrfj, is shown as a function of 
momentum for e-like and /i-like events. The value of the asymmetry for multi-GeV 
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Fig. 2. Angular distributions for e-like (left) and //-like (right) events, for sub-GeV (top) and multi- 
GeV (bottom) samples. The bars show the MC no-oscillation prediction with statistical errors, 
and the line shows the oscillation prediction for the best-fit parameters, sin 2 28 — 1.0 and Am 2 = 
3.5 x 10- 3 eV 2 . 
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Fig. 3. Angular distributions for e-like (left) and //-like (right) multi-GeV events (FC+PC), divided 
into 10 angular bins. 
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Fig. 4. Asymmetry as a function of momentum for e-like and /i-like events. For PC events the 
average neutrino energy is 15 GeV (the produced lepton exits and does not deposit all of its energy. 
The dashed line shows the oscillation prediction for the best-fit parameters. 



events (including PC) is A data = —0.31 ± 0.04, more than 7a from the expected 
value of zero, Amc = 0.01 ± 0.01. 



3.2.4. Oscillation Fit 

The oscillation fit is performed by constructing a x 2 quantity, 

X = L ~2 + ^^2" 

p,cos9,e,fj, j j 

The sum is over five bins in cos 9, seven bins in logp, where p is momentum, and two 
particle ID bins (e and p). Ndata is the number of data events in the bin, a is the 
statistical error on data and MC, and N osc is a weighted sum over events: 

£data 



N o->=T7Z £ Msin 2 2£,Am 2 ,a, ei ), (7) 
events 



C>mc 



where a is the flux normalization, sin 2 29 and Am 2 are the oscillation parameters, and 
£-data and Cmc are the data and MC livetimes. The ej values and their correspond- 
ing a/s take into account systematic uncertainties. Table |3| lists these systematic 
parameters, their estimated sigmas and their best-fit values. 
To build an allowed region: 



Table 3. Estimated systematic parameters. Note that the flux normalization a is not constrained 
in the fit. 



Systematic parameter 


Estimated a 


Value at best-fit 


Flux normalization, a 


(~ 25 %) 


8.4% 


Energy spectrum power law, £ ,2 - 7±<5 


0.05 


-0.004 


/i-like/e-like ratio, sub-GeV, j3 s 


8% 


1.5% 


/U-like/e-like ratio, multi-GeV, (3 m 


12% 


-12.5% 


FC/PC ratio, p 


8% 


-2.9% 


Zenith angle asymmetry, sub-GeV, i] s 


2.4% 


-0.4% 


Zenith angle asymmetry, multi-GeV, r/ m 


2.7% 


-0.2% 


< L/E u >, A 


15% 


-1.1% 



• This x 2 quantity is minimized at each sin 2 29, Am 2 value, with respect to the 
systematic parameters €j. 

• The minimum x 2 m sin 2 20, Am 2 parameter space, Xmin is found. 

• A confidence level map is built as a function of x 2 ~ Xmin- The unphysical 
region is taken into account according to a two-dimensional generalization of 
the Particle Data Book prescription OJQ). 

Figure |5] shows the allowed regions for 68%, 90% and 99% C.L. for — > u T 
oscillation. For this case, 90% C.L. corresponds to Xmm + 5-3. The best-fit oscillation 
parameter values in the physical region are sin 2 29 = 1.0 and Am 2 = 3.5 x 10~ 3 eV 2 , 
corresponding to x 2 /d-o.f. = 62/67. Figure |^ shows the distribution of residuals for 
the 70 terms in the x 2 sum at Xmin'i ^ is well fit by a normalized Gaussian. All of the 
best-fit €j values lie within their estimated uncertainties at the 1 a level. The x 2 /d.o.f. 
value for the no-oscillation hypothesis is 175/69, corresponding to a probability of less 
than 0.0001%. 

Figure [7] shows the new allowed region compared with the June 1998 region and 
the old Kamiokande resuK®. 

3.3. Upward-going Muons 

We can get neutrino oscillation information from the upward-going muon (upmu) 
sample by comparing the angular distribution to the expected one. In addition, 
although we can get no direct energy information, we do know that the stopping and 
through-going samples come from different parent neutrino energy distributions, so 
that the stopping to through-going ratio is a way of testing the energy dependence of 
the oscillation hypothesis. 

The through-going upward muon angular distribution is shown in Figure The 
overall stop/through ratio is given by = o.22 ± 0.03, compared to the expected 
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Fig. 5. Allowed regions for — ► v T oscillations, for FC and PC data. 
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Fig. 6. Fit residuals: distribution of \ 2 residuals for the 70 bins in the equation || sum at Xmin- The 



table on the left indicates the tj best-fit values. 
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Fig. 7. 90 % C.L. allowed regions compared with the June 1998 Super-K result and the Kamiokande 
result. 



value (using the Honda fluxB>) of = 0.37 ± 0.04 for the no-oscillation case. The 
stop/through ratio as a function of angle is shown in Figure [8[ 

The oscillation fit is done by construction of a x 2 similar to equation |^: 

2 ^ (®data ~ <$>oscY , f a \ 2 f T] ^ ' 

x =2^ — 



, i a- \(j a J \ J 



where the sum is over 10 zenith angle bins for through- going muons and five zenith 
angle bins for stopping muons. For this case, the weighted MC prediction is replaced 
by a flux calculation!* , and 

$ osc = $ ca £ C (sin 2 29, Am 2 ) * (1 + a) for through-going muons, and (9) 



&osc = < ^ ) caZc(sin 2 29, Am 2 ) * (1 + a)(l + 77) for stopping muons. (10) 

a is the absolute flux normalization (a a ~ 22%) and 77 corresponds to the relative 
stopping and through-going uncertainty (a v ~ 13%). The 90% C.L. allowed region 
(including both through-going and stopping upmu information) is shown in Figure |^. 
The best-fit oscillation parameters are sin 2 29 = 1 and Am 2 = 3.2 x 10~ 3 eV 2 , corre- 
sponding to x 2 /d.o.f. = 8/13. 

3.4- Combined Analysis 

We can combine FC/PC and upmu data© for a global fit by constructing a 
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Fig. 8. Left: Angular distribution of through-going upward muons. The circles represent the data, 
the solid line represents the normalized no-oscillation flux prediction, and the dashed line represents 
the best fit prediction for oscillations. Right: ratio of stopping to through-going muons as a function 
of cos 9. 
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Fig. 10. Allowed regions using combined information from FC, PC and upward-going muons. 



X 2 quantity including information from both the FC/PC and upmu x 2 quantities 
(equations | and ||): 

x 2 = T. {Ndata ~ 2 Nosc? +T.% (ii) 

where the sum is over the 70 FC/PC bins, and the 15 upmu bins. For FC/PC 
bins, N osc represents the weighted MC term of the FC/PC fit, whereas for the upmu 
bins, N osc represents the flux prediction. The systematic parameters €j include the 
same FC/PC terms as before; there is a common flux normalization parameter a 
(unconstrained, as for the FC/PC fit) and a common spectral index uncertainty S. 
In addition the parameter rj of the upmu fit is replaced by two parameters r]\ and rj2, 
each with an assigned a of 7%, representing the relative stopping and through-going 
uncertainty and the relative FC/PC and upmu uncertainty respectively. 

The resulting — > v T confidence level contours are shown in Figure 0. The 
best-fit value is at sin 2 29 = 1 and Am 2 = 3.2 x 1CT 3 eV 2 , corresponding to x 2 /d.o.f. 
of 70.2/82. 

3.5. FC Multiple Rings 



An additional FC sample can be used to check for consistency of oscillation behav- 
ior: the multiple ring component of the FC sample (all FC results so far have been for 



Zenith angle distribution of 2-ring, v e enriched sample 



Zenith angle distribution of 2-ring, v„ enriched sample 




Fig. 11. Left: angular distribution for v e enriched sample of multiple-ring events. The left-hatched 
bars show the MC prediction for no oscillation and the right-hatched bars show the prediction for 
the oscillation parameters sin 2 29 — 1 and Am 2 = 3.5 x 1CP 3 cV 2 . Right: ^-enriched sample. 

single rings, to tag the simplest quasielastic interactions). To produce enriched muon 
and electron flavor 2-ring samples, sub-GeV 2-ring events were required to satisfy the 
following conditions: 

• p > 300 MeV. 

• For z/ e -enriched: 2 e-like rings not in the ir° peak (see Section 4.2.2) or 1 e-like 
ring and 1 /i-like ring with the e-like ring having higher momentum. For either 
case, there must no muon decay electrons. 

• For z/^-enriched: 2 /i-like rings or 1 e-like ring and 1 /i-like ring with the /i-like 
ring having higher momentum. For either case, there must be at least one muon 
decay electron. 

According to MC, the enriched v e sample is 60% v e CC, 8% CC and 32% 
NC, and the enriched sample is 89% CC, 2% v e CC and 8% NC. The angular 
distributions for this preliminary analysis are shown in Figure [11]. Clearly, the results 
are consistent with the oscillation hypothesis. 

4. Oscillation Interpretation 

For all of the allowed regions so far, we have been assuming that the muon disap- 
pearance results from to v T oscillation. However, this assumption is not necessarily 
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Fig. 12. Angular distributions with the best fit for — > v e superimposed. 



true. What can we say about the flavors involved in the oscillation? 

4.1. z/ M v e 

Are the muon neutrinos oscillating into electron neutrinos? This is a fairly easy 
question to answer - the Super-K angular distribution by itself is inconsistent with 
pure — > v e oscillation. This assumption (including matter effects in the Earth) 
yields a very poor fit: x 2 /d-o.f. = 110/67 (see Figure |l|). In addition, the CHOOZ 
experiment in France rules out this hypothesis (for u e disappearance)!!^. 

4.2. -> v s 

Could the muon neutrinos be oscillating into a non-weakly interacting sterile neu- 
trino? With Super-K data alone, it is hard to tell this hypothesis apart from the 

— > v T hypothesis; we expect only tens of CC v T interactions in the current sam- 
ple, and the products of such interactions in our detector are nearly indistinguishable 
from other atmospheric neutrino events. However, an oscillation to a sterile neutrino 
could nevertheless leave subtle imprints on some properties of our data set. Here two 
preliminary studies are described. 



4.2.1. Angular Distribution 

Due to the fact that sterile neutrinos do not exchange Z bosons with nucleons in 



600 



400 



200 



Super-Kamiokande 45 kt-yr Pre II 




sub-GeV jx-like 



200 



100 



multi-GeV e-like 







multi-GeV u-like + PC 



sirf28=1 
Am 2 =4.5x10" 3 eV : 



-0.8 -0.4 0.4 0.8 -0.8 -0.4 0.4 0.8 

cos© cos© 



Fig. 13. Angular distributions with the best fit for — > v s superimposed. 



the earth, the angular distribution for — > v s will be distorted with respect to the 
— > u T oscillation case (see, for instance, Lipari and Lusignoli 0)). This effect is 

expected to be largest for high E and low Am 2 cases. The effect is, however, smeared 

over a distribution of neutrino energies. 

Figure [13] shows the fit for — > u s : a good fit (x 2 /d-o.f. = 64/67) is obtained. 

The corresponding allowed region is shown in Figure it is slightly smaller than 

the corresponding — > v T region, due to additional constraints from the matter 

effect in the earth. We have not yet included the high energy upward muon sample 

in this analysis. 

4.2.2. Single Neutral Pions 

A rather promising way to distinguish the tau and sterile oscillation cases is via 
NC neutrino interactions in the detector which produce single neutral pions. Sterile 
neutrinos have no NC interactions; therefore, fewer ir° events than expected for the 

— > v T case would be a signature of oscillation to a sterile neutrino. 

The signature of a it NC event in the detector is a pair of e-like rings (7-induced 
electromagnetic showers) from the decay of the pion, tt° — > 77. We can select a 
sample of such events by selecting events with exactly two e-like rings and no muon 
decay electrons. Figure [15] shows the reconstructed invariant mass distribution of this 
selected sample. The ir° peak is clearly visible. We further reduce and enrich the sam- 
ple by making an invariant mass cut of 90 MeV/c 2 < M 77 < 190 MeV/c 2 . According 
to MC, the events surviving this cut are an 84% pure sample of NC interactions. The 
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Fig. 14. Allowed region for — ► v s . 



angular distribution of these events is shown in Figure To 



We can also take the ratio of 7r°/e-like events and compare it to the expected value 
for MC. The FC single ring e-like sample is mostly CC interactions of non-oscillated 
VgS. This ratio of ratios should be nearly 1 if muon neutrinos are oscillating to tau 
neutrinos (because tau neutrinos interact via NC), and less than 1 if muon neutrinos 
are oscillating into sterile neutrinos. Table f| shows the numbers. 



Table 4. n° candidates. 





Data 


MC 


7r u -like 

7r°-like, background-subtracted 


279 
231.8 


253.6 
195.9 



The result is: 

(n°/e) DATA 



1.03 ± 0.06(data stat) ± 0.02(MC stat) ± 0.24(sys) (12) 



(7r /e) M c7 

If the background is subtracted, the result is: 

{n°/e) D ATA 



(7r°/e) 



MC 



1.11 ± 0.06(data stat) ± 0.02(MC stat) ± 0.26(sys) (13) 



Notice that the systematic contribution to the uncertainty exceeds the statisti- 
cal contribution by a large factor. The systematic uncertainty is in fact dominated 




Fig. 15. Invariant mass distribution for 7r°-like events selected according to the criteria given in the 
text. The mass cut is shown. 



by uncertainties in the ir° production cross-section in water. The K2K collabora- 
tion @) hopes to reduce this uncertainty by high statistics measurements of neutrino 
interactions in its near water detector. 



5. Summary 

Table |5] summarizes the status of fits of the new Super-K data to three possible 
two-flavor oscillation scenarios. We obtain good fits for — > v T and — > v s \ pure 
two-flavor — > u e is ruled out. More work is currently in progress. 
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Table 5. Summary of two- flavor oscillation fits for new Super-K results. 



No-oscillation fits 








X 2 /d.o.f. 


FC/PC 
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